Effect of magnetic fields on the spin evolution of non-polarized 
87 Rb Bose-Einstein condensates 
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The spin mixing dynamics of spin-1 Bose-Einstein condensates with zero magnetization and under 
an external magnetic field is investigated. The time-dependent solutions are obtained via a diago- 
nalization of the Hamiltonian, which has a simple form under the single mode approximation. The 
features of evolution are compared in detail with those with the field removed so as to emphasize 
the effect of the field, which can induce strong oscillation in population of atoms in spin component 
. A new mode of oscillation characterized by a high frequency and a low frequency, is found when 
the field is sufficiently strong. 
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I. INTRODUCTION 

The rich physics in the Bose-Einstein condensate of 
atoms with spin degrees of freedom has attracted consid- 
erable interests in the last decade. [J 0, H H i, i, 0, HI 
Experimentally, one has been able to prepare spinor con- 
densates of given number of atoms and given magneti- 
zation, and then to observe their evolutions. The time- 
dependent population of each spin component has been 
measured by the Stern-Gerlach splitting technique. Var- 
ious modes of oscillation of the populations have been 
found, they are sensitive to the initial states. [1, d, 
Coherent dynamics has been observed experimentally. 

egi lata 

When an external magnetic field is applied, the evo- 
lution will be affected. This is a way to control the 
evolution. Efforts along this line have been made. 
EE EE El El E! HI Most theoretical calculations on 
this problem are based on the mean-field theory. In this 
approach, the spin dynamics can be modeled as a non- 
rigid pendulum which displays a variety of periodic os- 
cillations in magnetic fields. [H HlJ The dependence of 
the oscillation period on the magnetic field has been con- 
firmed by experiments. [22j In fact, theoretical studies 
suggested that more complicate spin dynamic behavior 
would be possible. [H [H[ 

In a previous paper of us [24| , the spin evolution with- 
out a magnetic field has been studied. Where, instead 
of using the mean field theory, the time-dependent solu- 
tions of the Hamiltonian have been obtained under the 
single mode approximation (SMA). Thereby an analyt- 
ical formula governing the evolution has been derived. 
The present paper is a direct generalization of the previ- 
ous one, the aim is to clarify the effect of magnetic field 
on the spin evolution. For this purpose, the Hamilto- 
nian arising from the SMA has been exactly diagonal- 



ized. Thereby the time-dependent populations of spin-0 
component starting from various initial states can be ob- 
tained. We have focused on the condensates with zero 
magnetization. Since the solutions of our Hamiltonian 
are exact, both short-term and long-term behaviors can 
be predicted by the theory. In the follows the basic 
approach is given in Section II. The features of spin- 
evolution without a magnetic field is briefly reviewed in 
Section III. Those under a magnetic field are presented 
in Section IV. In Section V, the effects of atom number 
and the confinement potential are discussed. Finally in 
Section VI, a summary is given. 



II. HAMILTONIAN AND ITS SOLUTION 

Let Uij — (co + c^Fi • Fj)S(Ti — r^) be the spin- 
dependent interaction among N spin-1 atoms. Under 
the SMA, all the atoms have the same spatial normal- 
ized wave function 4>(r). Accordingly, the Hamiltonian 
for the spin evolution under a magnetic field B with di- 
rection lying along the Z-axis reads 



H = gS 2 



(i) 
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where g = ^C2 / dr|</>(r)| 4 . The second term arises from 
the quadratic Zeeman effect, where q = (E + + E_ — 
2E Q )/2 is the energy difference of the three Zeeman levels, 
which is related to B as q = fi 2 B B 2 /(h 2 E HF s), and E HF s 
is the hyperfine splitting. The linear Zeeman term does 
not affect the dynamics due to the conservation of the 
total magnetization M. 

Let a Fock-state be denoted as \N\,Nq,N-\), where 
is the number of atoms in spin component fi. Since 
N±i = (N - N ± M)/2, due to the conservation of N 
and M which are initially given, the Fock-state can be 
simply denoted as \Nq) with No ranged from V — M, 
N — M — 2, toOorl. They form a complete set of bases 
for the many particle spin space. The matrix elements of 



2 



H read 



(N \H\N ) = [gA (M,No) + q(N-N )}5 No , ^ 
+ffA f (M,JVb)^ 0l ah+2 



- ff ^_(M,JV )J Woi w5 _2 



(2) 



where 



A (M, AT ) = M 2 + N + No + 2NN - 2A 2 



III. EVOLUTION WITHOUT MAGNETIC 
FIELD 

In order to understand the effect of magnetic field B, 
we first review the evolution with B = 0. The evolution 
is described by the formula of [24[ , which is exact for the 
Hamiltonian {T]) with q = 0. In order to obtain numerical 
results from the formula, N = 400 and ui = 3000/ sec 
are firstly assumed. Then the effect of N and uj will be 
studied. The evolution has the following features: 



A+(M,N ) = [N (N ~l){N+M~N +2)(N~M-Na+2)} 1 ^ 



A_(M,N ) = [(A + l)(A +2)(A+Af-Ar )(A-Af-A )] 1 / 2 

After the diagonalization of the Hamiltonian, the 
eigenenergies Ej and eigenstates \<fij) can be obtained, 

and we have H\cj)j) = Ej\(/>j), \<j)j) = J2n c tv \ n o)- lt is 
emphasized that both Ej and <pj are exact for the Hamil- 
tonian because the space for the diagonalization is com- 
plete. 

Let the initial state be a Fock-state with given N, M, 
and AT |t = o = I, and is denoted as |7). Then the total 
spin-state at time t is 



(i) Pf is 

•K/(\g\ui) and symmetric with respect to t p /2. 



strictly periodic with a period u p 



h \i) = Y J z- iE > t,h \4>i)(m 



-iEjt/h\ 



(3) 



It is emphasized that (3) is an exact time-dependent 
solution of the Hamiltonian, every detail of the evolu- 
tion can be thereby extracted. We are interested in the 
t— dependence of the population of spin-0 component, 

A 

namely, (^(t)\N \^(t))/N = P/ f (i), which can be de- 
rived from (3) as 



Pf(t) 



Bf 



Of(t) 



where 



(4) 



(5) 



No 



Of(i) = 2 ^ cos[(^ - EjWhftcl ^ c? No c> No N /N 

j<3' N 

(6) 

Eq.((4]) implies an oscillation around the background Bj . 

We consider the condensate of 87 Rb atoms trapped 
by a harmonic potential ^muj 2 r 2 . In the present pa- 
per, hoj, mG and sec are used as units. For 87 Rb, g ~ 
-6.57 x 10~ 5 (w/A 3 ) 1 / 5 (evaluated under the Thomas- 
Fermi approximation), and q ~ 1447T x 10~ 6 B 2 /uj. Nu- 
merical results with discussions are given in the follows. 



tp 

is anti-symmetric with respect to t p /A. Therefore, the 
study is sufficient to be confined in the duration (0, 
tp/A). Under the Thomas- Fermi approximation, t p « 
1.521tt(A/w 2 ) 3 / 5 x 10 4 sec (for 87 Rb). For the given N 
and uj, the above duration is 29.24sec. 

(ii) For zero-polarized condensates [M = 0), all the Pj 
have nearly the same ba ckg round Bj « 1/2 disregarding 
I as shown by eq. © of [2J| and by FigOJi of the present 
paper. 

(iii) If N is even, I must be even so that M can be 
zero. When 1 = 0, starting from the initial state that 
half atoms are spin-up while the other half are spin-down, 
P[J goes up directly to Bj rj 1/2 without oscillation as 
shown by the black curve in [T£i. From this curve we 
know that half of the ^ / atoms become n = within 
O.Olip, it certainly implies the occurrence of strong spin- 
flips | + 4< — > + 0. Then it remains extremely steady 
in a very long time until t is close to t p /4, where a round 
of strong oscillation occurs suddenly. 

(iv) When / = 2, a round of oscillation emerges in the 
early stage (red curve). When I = 2k, k rounds of oscil- 
lation emerge (blue curve). In general, Pj oscillates in 
the early stage but suddenly becomes very steady (this 
was first found in [5j), then the oscillation suddenly re- 
covers. Therefore there are zones of oscillation (ZOO) 
and zones of steady evolution (ZOS). They appear al- 
ternately, namely, ZOO-ZOS-ZOO as shown by the dark 
cyan curve, and again. In general, when I < N/2, a 
larger I leads to a broader ZOO and accordingly a nar- 
rower ZOS. Incidentally, if N is odd, I is odd. When 
I = 2k + 1 < N/2, the ZOO contains ft + 1/2 rounds of 
oscillation. 

(v) When I = N/4, N/8 rounds of oscillation are 
contained in the ZOO (dark cyan curve), It turns out 
that, in this occasion, the duration of the ZOO is just 
ip/12, that of the ZOS is also t p /12. Therefore the two 
ZOO and the one ZOS of the dark cyan curve all have 
the same duration t p /12. Meanwhile, the average fre- 
quency of oscillation in the ZOO is equal to 3N/(2t p ) « 
3.14(7Vu; 3 ) 2 / 5 x 10~ 5 , which w 5.13/ sec for our parame- 
ters. It can be imagined that, when N is very large, the 
frequency would be very high and the ZOO would appear 
as a band with a width. 

(vi) When I is close to N/2, the duration of the ZOO 
becomes very long (close to t p /8), and the amplitude of 
oscillation becomes very small except in a small domain 
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close to tp/8 as shown by the green curve in QJ,. Ac- 
cordingly, the ZOS in between becomes very narrow and 
unsteady 

(vii) When / = N/2, the two previous ZOOs trans- 
form to ZOSs, while the previous ZOS becomes to a very 
narrow ZOO containing only one round of oscillation sit- 
uated at tp/8 as shown by the top curve in[T^,. 

FigQJi is for the cases with / < N/2. The cases with 
/ > N/2 are given in FigHK- It was found that P° and 
Pjv_/ are one-to-one roughly similar, however they are 
greatly different when t is close to or t p /4. The early 
stage of evolution is referred to Fig[2K and0^,. 



IV. EVOLUTION UNDER A MAGNETIC FIELD 

When B ^ 0, via the procedure of diagonalization, we 
obtain the following numerical results. 
(I) Weak field 

The evolutions of P° are given in the (b) and (c) panels 
of FigQ]to[H we found 

(i) The effect of B is very sensitive if / is close to or 
N. Otherwise is less affected. Note that of the 
lowest three curves in QJi ([3k,) are much higher (lower) 
than their I/N. The effect of B is to pull them down 
(push them up) so that are closer to I/N. 

(ii) P° and P% are extremely sensitive to B. This 
is shown by the black curves in [TJa and [3Jd. The black 
curve representing Pq in Fig. la is greatly pulled down 
by the field when B = 10 ([lb), and it becomes nearly 
a horizontal line close to zero when B = 30 ((TJ:). Thus, 
for Pq, the spin- flips which occurs strongly in the early 
stage is severely suppressed by the field. 

For Pjy with all the atoms in fj, = initially, we 
know from FigJ3]that 64% (74%) of the atoms would be 
changed to fi ^ at 0.008i p ( 0.005i p ) if B = (B = 30). 
Thus, instead of being suppressed, the strong spin-flips 
in the early stage are accelerated by B. Consequently, a 
strong oscillation is thereby induced. When B=10, the 
oscillation is stronger in the duration (£ p /8,£ p /4), refer to 
the black and red curves of[3b- When B—30, the oscilla- 
tion is extremely strong through out nearly all the time 
(black curve of[5J;), this is a noticeable point and we will 
return to it. 

(iii) When / is neither close to nor N, the effect of 
B would be to shorten the ZOO and enlarge the ZOS 
if I < N/2 (e.g., the duration of the ZOO of the dark 
cyan curve of lc is £ p /13.2 instead of t p /12 in[]J,), or 
reversely if I > N/2 (e.g., refer to the pink curves in 
FigEl). Another effect of B is to spoil the stability of the 
ZOS, and cause "irregular" oscillation in both ZOO and 
ZOS. 

(iv) It seems that a new mode of oscillation might be 
caused by B. This is hinted, say, by the wavy pink curves 
in [3j3, [3}:, and the blue curve in[TJ;. This point is further 
studied below. 

(v) The strict periodicity and the inherent symmetry 
exist no more if B =/= 0. 



(II) Strong field 

The early stage of evolution of Pj under a strong field 
is plotted in Fig[5] within the duration (0, f p /40). Com- 
paring 5a with Figf2fc, we found that the curves with I 
close to zero are further suppressed so that P^ ~ I/N to- 
gether with a negligible oscillation. Comparing 5a with 
FigHh, the curves with / close to N keep their oscillation. 
The curves with / neither very small nor very large are 
relatively less affected, and they become more or less sim- 
ilar to each other as shown by the blue, pink, and navy 
curves ofEK For these curves the division into ZOO and 
ZOS holds no more. Bj of all the curves in 5a are close 
to I/N. 

When the field is even stronger as in[5jD, a new mode 
of oscillation emerges, where the new patterns contain a 
series of pulses, each has an olivary shape and includes 
certain rounds of oscillation inside. The number of round 
would decrease if I increases. The amplitudes of the oscil- 
lation is larger if / w N/2, and it will be too small to be 
seen if / « or N. li B increases further, the olives would 
become longer and thinner, and include more rounds. In 
the strong B limit the oscillation is completely suppressed 
and all P° become just horizontal lines. 

When B is sufficiently large, the first term of the 
Hamiltonian gS 2 can be considered as a perturbation. 
The set of the first order perturbative solutions of the 
Hamiltonian is 

Wi) = \I) + ^ q A-{M,I)\I + 2)--^A + (M,I)\I-2) (7) 

With this set, ^f(t) (eq.®) can be calculated analytically. 
Accordingly, we have 

Pf(t) « I/N+^[A 2 _(M,I)(l-cosat) 

-A\{M,I){1 -cos a'*)] (8) 

where a = [2q-g(4N- 81- 6)}u, a' = a- 16gw. Eq.® 
and © holds if \±A±(M,I)\ « 1. 

When M = and J = 0, A + (0, 0) = and A_(0, 0) = 
V2N, thus 

Pg(i) «2^(l-cosai) (9) 

When M = and I = N, A-(0,N) = and 
A + (0,N) = 2N, thus 

P° r (t)Ml_4^f-(l- cosa't) (10) 
T 

Both (J9j) and (fT0|) imply a small oscillation around a 
horizontal line as shown in [SJa, where the amplitude of 
the black curve is too small to be seen. 

When I is neither close to nor N — M, the sum of 
A 2 _(M, I) and A 2 _(M, I) is much larger than their differ- 
ence. Thus we have 

P°(i) « i/ N+ J^[A 2 _(0,I)-Al(0,I)-(A 2 _(0,I) 
+A 2 + (0, 1)) sin(^±^t) sin(^^i)] (11) 
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Since a + a' is much larger than a — a' , (fTTj) describes 
a high frequency oscillation included in a low frequency 
oscillation as shown in Fig 03 ■ This explains the origin of 
the olivary shapes. The factor sin( a ~" t) has a period 
7r/(4|<?|aj). Therefore the length of the olive would tend 
to 7r/(8|<7|w) = tp/8. Furthermore, Pj(i) would oscillate 

around I/N+^{A 2 _(Q, I)-A\(0, I)}, where the second 
term is small, as shown in and in FigfSb. 

V. THE EFFECT OF THE CONFINEMENT 
POTENTIAL AND THE PARTICLE NUMBER 

The above numerical results are obtained with u> = 
3000 and N = 400. From the formula of our previous 
paper [24| we know that the change of ui is simply equiv- 
alent to a change of the scale of time, i.e., all the curves 
of Pj (t) would remain unchanged except that the impli- 
cation of t p is changed with u> via the relation t p oc cj~ 6 / 5 
. On the other hand, since g oc lu 1 ^ , the decrease of u> 
will slightly increase the ratio q/g, therefore will amplify 
slightly the effect of B. 

When N changes, the period is also changed due to 
t p oc N 3 / 5 . When the time is re-scaled, the effect of 
N on the evolution is mild if B = 0. This is shown in 
FigG)h where N = 4000. Fig[5^ is very similar to FigQJ, 
although N is ten times larger (note that the time scales 
of these two figures are different). In particular, the ZOO 
and ZOS exist, and the rounds of oscillation contained in 
the ZOO is also 1/2. Furthermore, the ZOO of the dark 
cyan curve in [5^ has also the duration (0,i p /12), how- 
ever 500 rounds of oscillations arc included inside. Since 
g ex iV~ 3 / 5 , the increase of N will increase the ratio q/g, 
therefore will also amplify the effect of B. Comparing^ 
with []J;, the black curve is more severely suppressed, the 
red and blue are more severely pulled down. 

VI. SUMMARY 

The effect of the magnetic field on the spin-evolution of 
zero-polarized 87 Rb condensates has been studied. The 
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following points are summarized: 

(i) Periodicity and inherent symmetry. This important 
feature of Pj 1 , is spoiled by the field. 

(ii) Background. The field pushes all the B° from 
w 1/2 towards I/N. 

(iii) Zone of steady evolution. The ZOS is no more 
highly steady. When B is stronger, the division into ZOO 
and ZOS does not hold. 

(iv) Sensitivity. When I is close to or N, P° is highly 
sensitive to B. Otherwise, it is less sensitive. This is fur- 
ther shown in Fig(7j From this figure one can see that B 
causes oscillation in general. The frequency of oscillation 
will increase if B increases. The amplitude will firstly in- 
crease with B. However, when B exceeds certain values, 
the amplitude decreases with B. When B is sufficiently 
large, the amplitude will become so small that P® looks 
just like the horizontal line I/N. This situation would 
occur if B > 60, 320, and 2000 mG, respectively, when 
I = 0,N, and N/2. It implies that the suppression of the 
amplitude is much easier to be fulfilled if I w or N. On 
the other hand, in[7J: with / = N/2, P%/ 2 ^ s om y weakly 
disturbed even if B is as large as 50mG (red curve) , and 
the amplitude is still not very small even if B is as large 
as 450mG (orange curve). 

(v) Strong oscillation. When / is close to N, a mediate 
field (roughly, 30 < B < 180 mG) causes a particularly 
strong oscillation of P" through out nearly all the time 
as clearly shown in FigJTJs . This is a noticeable point. 

(vi) New modes of oscillation. FigJTJi shows clearly 
how the increase of B leads to the appearance of the new 
mode with olivary shape of pulses. That is similar to the 
quantum beats. 
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FIG. 1: (colored online) Evolution of P?(t) with various I 
and B. All figures in this paper are for 87 Rb atoms with 
uj = 3000 and N = 400 (the only exception of N is Fig.6). 
Accordingly, the period t p = 117sec, and t is given from 
to tp/4 (the case with t > t p /4 can be understood from the 
inherent symmetry). I is from to N/2 marked by the curves. 
Each curve has been shifted up by 0.2 more than its lower 
neighbor. 
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FIG. 2: (colored online) The same as Fig[T]but given only in 
the early stage. 




FIG. 3: (colored online) The same as Fig[T] but with I from 
N/2 to N. 
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FIG. 4: (colored online) The same as Fig|3]but given only in 
the early stage. 




FIG. 5: (colored online) Early stage of evolution of P?(t) 
under a strong B. I is from to JV marked by the curves. 
Each curve has been shifted up by 0.1 more than its lower 
neighbor, t is from to t p /40. 
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FIG. 6: (colored online) Similar to Figfflbut with N = 4000. 
t is from to t p /8. A new set of / each is marked by the 
associated curve. 
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FIG. 7: Evolution of P?(t) with various / and B. B is marked 
by the curves. The dotted lines are corresponding to B = 
. Each curve of (b) and (c) has been shifted up by 0.1 more 
than its lower neighbor, t is from to i p /4 (a and b) or to 
tp/40 (c). 



